Magnetic tight-binding and the iron-chromium enthalpy anomaly

We describe a self consistent magnetic tight-binding theory based in an expansion of the Hohenberg-Kohn density functional to second order, about a non spin polarised reference density. We show how a first order expansion about a density having a trial input magnetic moment leads to the Stoner--Slater rigid band model. We employ a simple set of tight-binding parameters that accurately describes electronic structure and energetics, and show these to be transferable between first row transition metals and their alloys. We make a number of calculations of the electronic structure of dilute Cr impurities in Fe which we compare with results using the local spin density approximation. The rigid band model provides a powerful means for interpreting complex magnetic configurations in alloys; using this approach we are able to advance a simple and readily understood explanation for the observed anomaly in the enthalpy of mixing.Comment: Submitted to Phys Rev

Forcing boundary-layer transition on an inverted airfoil in ground effect

The influence of the laminar boundary-layer state on a wing operating in ground effect has been investigated using experiments with a model that provides two-dimensional flow. The effect of a boundary-layer trip placed at varying distances from the leading edge was observed at various incidences in terms of on-surface characteristics, including pressure measurements, flow visualization, and hot-film anemometry, and off-surface characteristics with velocity surveys below and behind the wing. The act of forcing transition led to downforce being reduced and drag increased, moreover, it altered almost all aspects of the wing’s aerodynamic characteristics, with the effect becoming greater as the trip was placed closer to the leading edge. These aspects include the replacement of a laminar separation bubble with trailing-edge separation, a thicker boundary layer, and a thicker wake with greater velocity deficit. The importance of considering laminar phenomena for wings operating in ground effect has been shown

Glycaemic control targets after traumatic brain injury: a systematic review and meta-analysis.

BACKGROUND: Optimal glycaemic targets in traumatic brain injury (TBI) remain unclear. We performed a systematic review and meta-analysis of randomised controlled trials (RCTs) comparing intensive with conventional glycaemic control in TBI requiring admission to an intensive care unit (ICU). METHODS: We systematically searched MEDLINE, EMBASE and the Cochrane Central Register of Controlled Trials to November 2016. Outcomes of interest included ICU and in-hospital mortality, poor neurological outcome, the incidence of hypoglycaemia and infective complications. Data were analysed by pairwise random effects models with secondary analysis of differing levels of conventional glycaemic control. RESULTS: Ten RCTs, involving 1066 TBI patients were included. Three studies were conducted exclusively in a TBI population, whereas in seven trials, the TBI population was a sub-cohort of a mixed neurocritical or general ICU population. Glycaemic targets with intensive control ranged from 4.4 to 6.7 mmol/L, while conventional targets aimed to keep glucose levels below thresholds of 8.4-12 mmol/L. Conventional versus intensive control showed no association with ICU or hospital mortality (relative risk (RR) (95% CI) 0.93 (0.68-1.27), P = 0.64 and 1.07 (0.84-1.36), P = 0.62, respectively). The risk of a poor neurological outcome was higher with conventional control (RR (95% CI) = 1.10 (1.001-1.24), P = 0.047). However, severe hypoglycaemia occurred less frequently with conventional control (RR (95% CI) = 0.22 (0.09-0.52), P = 0.001). CONCLUSIONS: This meta-analysis of intensive glycaemic control shows no association with reduced mortality in TBI. Intensive glucose control showed a borderline significant reduction in the risk of poor neurological outcome, but markedly increased the risk of hypoglycaemia. These contradictory findings should motivate further research

A simple environment-dependent overlap potential and Cauchy violation in solid argon

We develop an analytic and environment-dependent interatomic potential for the overlap repulsion in solid argon, based on an approximate treatment of the non-orthogonal Tight-Binding theory for the closed-shell systems. The present model can well reproduce the observed elastic properties of solid argon including Cauchy violation at high pressures, yet very simple. A useful and novel analysis is given to show how the elastic properties are related to the environment-dependence incorporated into a generic pairwise potential. The present study has a close link to the broad field of computational materials science, in which the inclusion of environment dependence in short-ranged repulsive part of a potential model is sometimes crucial in predicting the elastic properties correctly.Comment: 10 pages, 3 figure

The stabilizing role of itinerant ferromagnetism in inter-granular cohesion in iron

We present a simple, general energy functional for ferromagnetic materials based upon a local spin density extension to the Stoner theory of itinerant ferromagnetism. The functional reproduces well available ab initio results and experimental interfacial energies for grain boundaries in iron. The model shows that inter-granular cohesion along symmetric tilt boundaries in iron is dependent upon strong magnetic structure at the interface, illuminates the mechanisms underlying this structure, and provides a simple explanation for relaxation of the atomic structure at these boundaries.Comment: In review at Phys. Rev. Lett. Submitted 23 September 1997; revised 16 March 199

An Effective-Medium Tight-Binding Model for Silicon

A new method for calculating the total energy of Si systems is presented. The method is based on the effective-medium theory concept of a reference system. Instead of calculating the energy of an atom in the system of interest a reference system is introduced where the local surroundings are similar. The energy of the reference system can be calculated selfconsistently once and for all while the energy difference to the reference system can be obtained approximately. We propose to calculate it using the tight-binding LMTO scheme with the Atomic-Sphere Approximation(ASA) for the potential, and by using the ASA with charge-conserving spheres we are able to treat open system without introducing empty spheres. All steps in the calculational method is {\em ab initio} in the sense that all quantities entering are calculated from first principles without any fitting to experiment. A complete and detailed description of the method is given together with test calculations of the energies of phonons, elastic constants, different structures, surfaces and surface reconstructions. We compare the results to calculations using an empirical tight-binding scheme.Comment: 26 pages (11 uuencoded Postscript figures appended), LaTeX, CAMP-090594-

Two-band second moment model and an interatomic potential for caesium

A semi-empirical formalism is presented for deriving interatomic potentials for materials such as caesium or cerium which exhibit volume collapse phase transitions. It is based on the Finnis-Sinclair second moment tight binding approach, but incorporates two independent bands on each atom. The potential is cast in a form suitable for large-scale molecular dynamics, the computational cost being the evaluation of short ranged pair potentials. Parameters for a model potential for caesium are derived and tested

The MOLDY short-range molecular dynamics package

We describe a parallelised version of the MOLDY molecular dynamics program. This Fortran code is aimed at systems which may be described by short-range potentials and specifically those which may be addressed with the embedded atom method. This includes a wide range of transition metals and alloys. MOLDY provides a range of options in terms of the molecular dynamics ensemble used and the boundary conditions which may be applied. A number of standard potentials are provided, and the modular structure of the code allows new potentials to be added easily. The code is parallelised using OpenMP and can therefore be run on shared memory systems, including modern multicore processors. Particular attention is paid to the updates required in the main force loop, where synchronisation is often required in OpenMP implementations of molecular dynamics. We examine the performance of the parallel code in detail and give some examples of applications to realistic problems, including the dynamic compression of copper and carbon migration in an iron-carbon alloy

Angular Forces Around Transition Metals in Biomolecules

Quantum-mechanical analysis based on an exact sum rule is used to extract an semiclassical angle-dependent energy function for transition metal ions in biomolecules. The angular dependence is simple but different from existing classical potentials. Comparison of predicted energies with a computer-generated database shows that the semiclassical energy function is remarkably accurate, and that its angular dependence is optimal.Comment: Tex file plus 4 postscript figure

Development of a tight-binding potential for bcc-Zr. Application to the study of vibrational properties

We present a tight-binding potential based on the moment expansion of the density of states, which includes up to the fifth moment. The potential is fitted to bcc and hcp Zr and it is applied to the computation of vibrational properties of bcc-Zr. In particular, we compute the isothermal elastic constants in the temperature range 1200K < T < 2000K by means of standard Monte Carlo simulation techniques. The agreement with experimental results is satisfactory, especially in the case of the stability of the lattice with respect to the shear associated with C'. However, the temperature decrease of the Cauchy pressure is not reproduced. The T=0K phonon frequencies of bcc-Zr are also computed. The potential predicts several instabilities of the bcc structure, and a crossing of the longitudinal and transverse modes in the (001) direction. This is in agreement with recent ab initio calculations in Sc, Ti, Hf, and La.Comment: 14 pages, 6 tables, 4 figures, revtex; the kinetic term of the isothermal elastic constants has been corrected (Eq. (4.1), Table VI and Figure 4